Measurement of particle and bubble accelerations in turbulence 
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We use an extended laser Doppler technique to track optically the velocity of individual particles 
in a high Reynolds number turbulent flow. The particle sizes are of the order of the Kolmogorov scale 
and the time resolution, 30 microseconds, resolves the fastest scales of the fluid motion. Particles 
are tracked for mean durations of the order of 10 Kolmogorov time scales. The fastest scales of the 
particle motion are resolved and the particle acceleration is measured. For neutrally buoyant par- 
ticles, our measurement matches the performance of the silicon strip detector technique introduced 
at Cornell University 0, This reference dynamics is then compared to that of slightly heavier 
solid particles (density 1.4) and to air bubbles. We observe that the acceleration variance strongly 
depends on the particle density: bubbles experience higher accelerations than fluid particles, while 
heavier particles have lower accelerations. We find that the probability distribution functions of 
accelerations normalized to the variance are very close although the air bubbles have a much faster 
dynamics. 

PACS numbers: 47.27. Jv,47.27.Gs,02.50.-r 



The Lagrangian approach to fluid dynamics is a nat- 
ural one when one addresses problems related to mixing 
and transport [ij. It has also been widely studied in the 
context of intermittency in fully developed turbulence. In 
recent years, several novel experimental techniques have 
been developed. The pioneering optical tracking method 
developed in the Cornell group has revealed that fluid 
particles experience extremely intense accelerations, with 
probability density functions (PDFs) having stretched 
exponential tails 0, Initially limited to very short 
particle tracks, the technique has been extended with 
the use of ultrafast optical cameras and is currently 
applied to the study of multiple particle statistics Q . In- 
dividual particles have been tracked for time duration of 
the order of the flow integral time scale using an acoustic 
technique Q: in an insonified volume, individual parti- 
cles scatter a sound wave whose Doppler shift carries the 
tracer velocity. In reason of the very fast decrease of the 
acoustic scattering cross-section, this method is limited 
to particles with size of the order of the wavelength, i.e. 
inertial ranges sizes when using acoustics. However 
the principle of the technique is completely analogous 
to laser Doppler velocimetry (LDV), provided that ex- 
panded light beams are used. Interference fringes are 
created at the intersection of two wide laser beams; a 
particle that crosses these fringes scatters light with a 
modulation frequency proportional to its velocity com- 
ponent perpendicular to the direction of the fringes Q. 
The advantage, compared to the acoustic method, is that 
the much smaller wavelength allows a better resolution 
in space and also the use of smaller tracer particles. 

In this Letter, we describe the principles of this tech- 
nique and validate it against the known features of the 
Lagrangian acceleration statistics in a fully turbulent 
von Karman flow at R\ < 850. We then apply it to 



track the dynamics of particles whose density differs from 
that of the fluid. The dynamics of such inertial par- 
ticles is relevant for many engineering applications re- 
lated to transport, mixing, dispersion, etc 0]. Signifi- 
cant theoretical and numerical progress in this domain 
has been made in the limiting case of infinitely heavy, 
pointwise particles 11| and has received experimental 
verifications Bin. We report the first experimental 
measurements of accelerations of particles having a den- 
sity in the range 10"^ (air bubbles) to 1.4 (PMMA) in 
the same highly turbulent flow. Taking into account the 
added mass effect for small spherical particles (i.e. the 
displacement of fluid elements by the particle motion), 
the effective density of the bubbles is only 3 times less 
than that of the fluid, while the PMMA particles are 
roughly 1.5 times larger. Because of the mismatch of 
density, light particle tend to be trapped in high vortical 
regions: as result of a lower inertia, the centrifugal force 
can not compensate the pressure gradient which drive 
them into the core of the vortices. On the contrary the 
centrifugal force is stronger than the pressure gradient 
for heavy particles so that they are ejected form vortex 
cores and concentrate in high strain regions [l3l ■ Because 
of this distinct spatial sampling of the flow, particle with 
different buoyancy are expected to exhibit quite differ- 
ent dynamical behavior. Indeed, we do find that the 
particles have different dynamical characteristics such as 
acceleration variance or correlation time. The PDFs of 
their accelerations remain close for value less than about 
10 times the acceleration variance, and differ for higher 
values. 

The Laser Doppler technique is based on the same prin- 
ciple as the ultrasound Doppler method which was shown 
to be very valuable for Lagrangian particle tracking Q . I 
order to access dissipative scales, and in particular for ac- 
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FIG. 1: Experimental setup, (left): schematics of the von 
Karman flow in water - side view, (right): principle of the 
Laser-Doppler Velocimetry using wide beams (ELDV) - top 
view of the experiment. PM: location of the photmultipler 
which detect scattering light modulation as a particle crosses 
the interference pattern created at the intersection of the laser 
beams. 



celeration measurements, one adapts the technique from 
ultrasound, to Laser: the gain is of a factor 1000 in wave- 
length so that one expects to detect micron-sized parti- 
cles. For a Lagrangian measurement, one has to be able 
to follow the particle motion to get information about its 
dynamics in time. For this, wide Laser beams are needed 
to illuminate the particle on a significant fraction of its 
path. The optical setup is an extension of the well known 
laser Doppler velocimetry technique; cf. Fig.l. A Laser 
beam is split into two beams; each is then expanded by 
a telescope so that their diameter is about 5mm. Then 
the two beams are directed in the flow. In their intersec- 
tion volume, they create an array of interference fringes. 
As a particle goes across the fringes, the scattered light 
is modulated and the frequency modulation is directly 
proportional to the component of the velocity perpen- 
dicular to the interference fringes. One then measures 
one component of the particle velocity. In practice, we 
use a CW YAG laser of wave length 532 nm with 1.2 W 
maximum output power. In order to get the sign of the 
velocity, the standard method consists in using a acousto- 
optic modulator (AOM) to shift the frequency of one of 
the beams so that the fringes are actually travelling at 
a constant speed. Here we use one AOM for each beam, 
the two excitation frequencies of the AOM being shifted 
by 100 kHz. The angle of the two beams is tuned to im- 
pose a 60 microns interfringe so that lOOkHz corresponds 
to 6 m/s. As the beams are not focused, the interfringe 
remains constant across the measurement volume whose 
size is about 5 x 5 x 10 mm^. The measurement volume 



FIG. 2: Signal processing, (top): raw light modulation, 
as detected by the photodiode when a particle crosses the 
fringes, (middle) : heterodyne detection of the frequency mod- 
ulation, (bottom): velocity trace extracted using the AML 
algorithm [Tsj . 



is imaged on a photomultiplier whose output is recorded 
using a National Instrument PXI-NI5621 digitizer. 

The flow is of the Von Karman kind as for the ultra- 
sound measurements Water fills a cylindrical con- 
tainer of internal diameter 15 cm, length 20 cm. It is 
driven by two disks of diameter 10 cm, fitted with blades. 
The rotation rate is fixed at values up to 10 Hz. For the 
measurements reported here, the Taylor based Reynolds 
number up to 850 and the dissipation rate e up to 25 
W/kg. We study three types of particles: neutrally buoy- 
ant polystyrene particles of size 31 microns and density 
1.06, PMMA particles of size 43 microns and density 1.4 
and air bubbles with a size of about 100 microns. The 
size of the bubbles, measured optically by imaging the 
measurement volume on a CCD, is imposed by the atom- 
isation of a large bubble by the turbulent flow, a process 
known to lead to a well defined and stationary size dis- 
tribution. 

The signal processing step is crucial as both time and 
frequency - (i.e. velocity - resolutions rely on its per- 
formance. Frequency demodulation is achieved using the 
same algorithm as in the acoustic Doppler technique. It 
is a approximated maximum likelihood method coupled 
which a Kalman filter [l^: a parametric estimator as- 
sumes that the signal is made of a modulated complex 
exponential and Gaussian noise. The amplitude of the 
sound and the modulation frequencies are assumed to be 
slowly evolving compared to the duration of the time win- 
dow used to estimate the instantaneous frequency. Here 
the time window is about 30 fis long and sets the time 
resolution of the algorithm. Outputs of the algorithm are 
the instantaneous frequency, the amplitude of the mod- 
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FIG. 3: (a) Probability density functions (PDFs) of the ac- 
celeration, normalized by its variance, for neutral particles, 
(symbols): ELDV measurements; (black thick line): Cornell 
data at R\ = 690. (b) autocorrelation coefficients of the ac- 
celerations for the same cases. 

ulation and a confidence estimate which can be used to 
discriminate bad data. An example of the light scattered 
by a particle is displayed in Fig. 2. together with the 
Dopplcr frequency modulation, and final velocity signal. 
Afterwards, the acceleration of the particle is computed 
by differentiation of the velocity output, using a low pass 
filter to smooth the noise as in Q. After processing, the 
data consists in a collection of sequences whose lengths 
are exponentially distributed. 

We first compare the data for neutral particles to the 
existing data recorded in a similar flow in Cornell Uni- 
versity using high speed imaging 0, Q : the measurement 
used linear cameras of 512 pixels running at speeds up 
to 70, 000 frames per seconds to get the time and lin- 
ear resolution required to resolve the acceleration. The 
probability density functions (PDFs) of the acceleration, 
from our ELDV technique at increasing Reynolds num- 
bers and for the Cornell data at R\ = 690 are displayed in 
Fig. 3. The distributions are seen to be very good agree- 
ment both qualitatively and quantitatively. In order to 
compare the two experiments, the following Heisenberg- 
Yaglom scaling is used 

{a')=aoe'/':.-'/\ (1) 

where e is the energy dissipation rate per unit mass and 
1/ ~ 1.3 • 10~^ m^.s~^ is the kinematic viscosity of the 
fluid. We derive here qq = 6.4 ± 1 at i?A = 850 compared 
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TABLE I: Parameters of the flow. SI: rotation rate of 
the disks, e dissipation rate obtained from the power con- 
sumption of the motors (with an accuracy of about 20%). 
The Tay lor-based tu rbulent Reynolds number is computed as 
Rx = ^/ibufrns/^i^! B.nd ao is derived from the Heisenberg- 
Yaglom relationship - Eqn. 1. 



to 6.2 ± 0.4 for the CorneU data at Rx = 690. The accel- 
eration variance is computed using the same procedure 
as in 3 : it is obtained for several width of the smooth- 
ing kernel used in the diiferentiation of the velocity signal 
and then interpolated to zero filter width. The accelera- 
tion autocorrelation function, shown in Fig. 3b, decays in 
a time of the order of the Kolmogorov time = y/iy/e, 
the fastest timescale of the turbulent flow - for the data 
shown here, = 0.25ms. This confirms that our tech- 
niques achieves a fast enough time and velocity resolution 
to get a good estimate of the acceleration. 

The very close agreement between the measurements 
reported here and the Cornell 0] data is expected be- 
cause they share the same flow geometry. The fact that 
it is observed using two very different techniques and sig- 
nal processing validates both (as the performance of the 
high speed imaging method had not been matched pre- 
viously). 

We now apply our technique to compare Lagrangian 
tracers to the dynamics of heavier and lighter particles. 
We flrst compute the velocity root mean square value 
Urms for all three cases: the values are 1.1,1.2, 1.0 ± 
O.lm.s-i at Rx = 850 for the neutral, heavy (PMMA 
spheres) and light particles (bubbles). Within error 
bars, the large scale dynamics seems to be unaffected 
by changes in the particle inertia. The acceleration dis- 
tribution and autocorrelation in the three cases is shown 
in Fig. 4. The acceleration PDFs are quite similar for 
moderate acceleration values (below about lOorms), as 
also observed in low Reynolds number numerical simu- 
lations jl^l. However, the probability of very large ac- 
celerations is reduced in the case of bubbles compared 
to neutral particles. The normalized acceleration vari- 
ance ao varies very significantly: it is reduced to 4.3 ± 1 
for heavier particles at Rx = 850 while it is increased to 
26 ± 5 for bubbles. The correlation functions also show 
significant changes with the inertia: the characteristic 
time of decay is longer for heavy particles and shorter 
for bubbles compared to neutral particles. We measure 
Tcorr/Tri = 0.5, 0.9, 0.25 respectively for neutral, heavy 
and light particles, with the correlation time defined as 
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FIG. 4; (a) Normalized autocorrelation coefficient of the ac- 
celeration, (b) Probability distribution function of accelera- 
tions, normalized to the variance of the data sets. Flow at 
i?A = 850. 



ties in the fluid 17[ but the boundary condition may be 
different from that of a solid particle. We find that the 
trend for the acceleration variance follows qualitatively 
that of (3 in eqn Quantitatively, for heavy particles 
it changes roughly as /3, but for bubbles, it is only about 
2.6 times that of the fluid. This is different from what 
was reported in [ist where a stronger influence of small 
change in the particles density was observed. 

To conclude, we have reported here an extended Laser 
Dopplcr velocimetry technique (ELDV) for the tracking 
of individual particles in fully developed turbulence. The 
advantage of our technique is that it may be more eas- 
ily adapted from commercial equipment than ultrafast 
PIV, and requires less laser power for illumination. Its 
application to the study of inertial particles with a flnite 
density with respect to the fluid shows that quite different 
dynamics may lead to very similar statistics of particle 
accelerations, an observation that may prove useful for 
modeling. 
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the half-width at mid amplitude of the correlation func- 
tion. We thus observe significant changes in the dynam- 
ics, even if the distribution of acceleration weakly changes 
with inertia. 

We now briefly discuss these results. In order to char- 
acterize the dynamics of a solid particle, one must specify 
two dimensionless numbers: the Stokes number which, in 
the case of turbulent flows, is the ratio of the response 
time of the particle to the Kolmogorov time scale, and 
the density ratio of the particle's material to that of the 
fluid. In the asymptotic case of very large inertia, only 
the Stokes number matters [ll| . Here the inertia remains 
flnite even for the bubble case because of the added mass 
contribution. In the limit of very small particles com- 
pared to the Kolmogorov length scale one can derive an 
effective equation of motion for a solid particle [lit : 



dvp 
It 



f3- 



Dt 



(2) 



where Vp is the particle velocity, /3 = ipf /{2pp + pf), pf 
and pp are the fluid and particle specific mass respec- 
tively), Du/Dt is the acceleration of the fluid particle 
that would be at the position of the solid particle in the 
undisturbed flow and F incorporates other forces such as 
drag, lift, history and possibly the buoyancy (negligible 
here). We have f3 = 1, 0.8 and 3 for the neutral, heavy 
and light particles. Bubble with such a small diameter 
are usually considered as being rigid because of impuri- 
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